Non-Abelian optical lattices: Anomalous quantum Hall effect and Dirac fermions 
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We study the properties of an ultracold Fermi gas loaded in an optical square lattice and subjected 
to an external and classical non-Abelian gauge field. We show that this system can be exploited as 
an optical analogue of relativistic quantum electrodynamics, offering a remarkable route to access 
the exotic properties of massless Dirac fermions with cold atoms experiments. In particular we show 
that the underlying Minkowski space-time can also be modified, reaching anisotropic regimes where 
a remarkable anomalous quantum Hall effect and a squeezed Landau vacuum could be observed. 
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Low energy excitations of fermionic lattice systems 
are usually governed by the non-relativistic Schrodinger 
equation. However, this description must be profoundly 
altered in the vicinity of Dirac points, where the en- 
ergy bands display conical singularities and quasiparti- 
cles become massless relativistic fermions. Such a re- 
markable behavior can be induced by a honeycomb ge- 
ometry [U m |3l [4[ |5], or by additional uniform [6 or 
staggered [71 [8] magnetic fields. Here we show that the 
natural playground for emerging Dirac fermions is pro- 
vided by multi-component fermionic atoms subjected to 
artificial non-Abelian gauge fields. We emphasize that 
these external fields can be produced by generalizing the 
recent experiment 0, as proposed in [TOl |TT]. Such 
gauge fields give rise to intriguing phenomena such as the 
non-Abelian Aharonov-Bohm effect fTO^, generation of 
magnetic monopoles [12 , non-Abelian atom optics [13], 
quasi-relativist ic effects [M], or even the modification of 
the metal-insulator transition [TS . In this Letter, we 
show that the physical properties of massless relativis- 
tic fermions are completely characterized by the non- 
Abelian features of the external gauge fields. Further- 
more the anisotropy of the underlying Minkowski space- 
time can be controlled externally, producing an anoma- 
lous quantum Hall effect characterized by a squeezed 
Landau vacuum. 

We consider a system of two-component (two-color) 
fermionic atoms trapped in an optical square lattice with 
sites at r = (n, m)a, where a is the lattice spacing and 
n^m ^Ij. In the non-interacting limit, which can be ob- 
tained by means of Feshbach resonances [T6^, fermions 
freely hop between neighboring sites. The addition of an 
external gauge potential A modifies the hopping Hamil- 
tonian according to the Peierls substitution 

= E E 4' (r') e-' ^-^^ c.(r) + h.c, (1) 

(r,r') tt' 

where t is the hopping amplitude, Cr{v) is the fermionic 




Figure 1: (a) Square lattice subjected to a non-Abelian 
gauge potential. This external field induces state-dependent 
hoppings described by the U(2) operators Ux and Uy. (b) 
Energy bands close to the vr-flux regime = 7r/2 + 0.1, 
— 71-/2 — 0.1), with vanishing Abehan flux $ = 0. The 
bands touch at four Dirac points inside the first Brillouin zone 
(BZ), where the energy scales linearly with momenta E ^ k. 



field operator in color component r = 1,2, and we set 
h = e = 1. Our setup features an external gauge 
potential with both commutative and non-commutative 
components A = ^{—y^x) + a{BaCry^ Bf^cix)^ where 
Bq^ 5q,, Bf^ are controllable parameters and (Jx^y are Pauli 
matrices. Accordingly the hoppings are accompanied 
by non-trivial unitary operators, Ux{m) = e~^^*^^e^*^"^^ 
and Uy{n) = e'^'^^e'^^''- , where ^ = Boa'^ is the Abelian 
magnetic flux and ^a,f3 = Bo^^^a^ are the non-Abelian 
fluxes (see Fig. [l]a). 



Let us point out that the gauge fields considered in this 
work can be realized following the proposals [lOl [TTl [17] , 
along the lines of the recent experiment [9^ , and provide 
non-Abelian analogues of homogeneous magnetic fields 
since they are characterized by constant Wilson loops. 
Indeed, atoms hopping around an elementary plaque- 
tte undergo a unitary transformation U = Ux{m)Uy{n -\- 
l)/7j(m + l)/7^(n), explicitly given by 



U 



^i27r<l> 



(cil + C2crz + csay + c^ax 



(2) 
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where the constants {cj} are hsted in [18]. For specific 
values of ^a,/3 the loop matrix reduces to a phase fac- 
tor and reproduces the Abelian 7r-flux = = f ) 
or Hofstadter {^^ = *^/3 = 0) models jSl |T9]. How- 
ever in general cases, it is a non-trivial U(2) opera- 
tor exhibiting non- Abelian properties such as the non- 
Abelian Aharonov-Bohm effect. The gauge- invariant 
Wilson loop W = tiU provides a clear distinction be- 
tween the Abelian (|W| =2) and non- Abelian {\W\ < 2) 
regimes. We stress that the Wilson loop is homogeneous 
and that the corresponding spectrum exhibits well devel- 
oped gaps [17]. 

In order to isolate non- Abelian effects, we first study 
the regime of vanishing Abelian flux ^ = 0. The 
Hamiltonian is diagonalized in momentum space and 
the fermion gas becomes a collection of non-interacting 
quasi-particles with energies shown in Fig{l|b. Close 
to the marginally Abelian regime {^a^^fd ^ '7r/2), the 
spectrum develops four independent conical singularities 
kD e {(0,0),(f,0),(0,f),(f,f)} G BZ, which corre- 
spond to massless relativistic excitations at half filling. 
Around these points p = k — ko, the low-energy proper- 
ties are accurately described by a Dirac Hamiltonian 

^eff = ^ ^j[,i^D^p, Hb = c^a^Px + CyayPy, (3) 
p 

where = (cip, C2p)^ is the relativistic spinor, the Dirac 
matrices a^^ay fulfill {aj^ak} = 26 jk (e.g. around ko = 
(0,7r/a), ax = cry and ay = dx)^ and Cx = 2atsin<l>Q,, 
Cy = 2atsm^i3 represent the effective speed of light. 
We stress here that the control over the non-Abelian 
fluxes ^a,i3 offers the exotic opportunity to modify the 
structure of the underlying Minkowski space-time, reach- 
ing anisotropic situations where Cx 7^ Cy. Hence, non- 
Abelian optical lattices provide a quantum optical ana- 
logue of relativistic QED, where the emerging fermions 
and the properties of the corresponding space-time rely 
on the non-Abelian features of the external fields. Fur- 
thermore, it is also possible to observe a transition be- 
tween relativistic and non-relativistic dispersion relations 
as the energy is increased. This abrupt change of the 
quasi-particle nature is revealed by Van Hove singular- 
ities (VHS) in the density of states, as displayed in 
Figs j2j a-c. 

The transport properties of 2D Fermi gases subjected 
to external gauge fields are characterized by the optical- 
lattice analogue of the well-known quantum Hall effect 
(QHE) [20 . In this context, the transverse Hall conduc- 
tivity measures the response of the system to a static 
force, e.g. a lattice acceleration, and takes on quantized 
values axy = ^ with G Z, when the Fermi energy Ep 
lies in a gap [17] . Surprinsingly, the quantized conductiv- 
ity of cold gases can be directly observed through density 
measurements thanks to the Streda formula f2T . Here 
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Figure 2: (a) Density of states (DOS) in the vr-flux regime 
<I>cK — $/3=7r/2 when $ = 0. (b) Hall conductivity in units of 
as a function of the Fermi energy in the same regime for 
$ = 1/41. Black arrows designate the VHS. (c) DOS close 
to the TT-flux regime = 7r/2 + 0.1 and = 7r/2 - 0.1 
when $ = 0. (d) Hall conductivity haxy — haxy^Ep) in the 
same regime as (c) for $ = 1/41. Dark red and light green 
arrows respectively designate the VHS E^^^ and ^^Je^ (cf. 
Eq. ([6|). 

we show that non-Abelian effects have dramatic conse- 
quences on the QHE which occurs when an additional 
Abelian flux <I> is applied to our system. The quantized 
values of the transverse conductivity are calculated as the 
sum of topological invariants associated to each energy 
band, the so-called Chern numbers [22j, 

= - E ^ / tr^(V;n)d/c, (4) 

where J^{i)n) = {dk^i^n\dkyi^n) - {dkyi^n\dk^i^n) is the 
Berry's curvature of the band En. Here the Chern num- 
bers are computed numerically by discretizing the Bril- 
louin zone [23 . A lattice gauge theory method allows to 
determine the Berry's curvature 

J'xyiki) = lnTx{ki)Ty{ki + x)Tx{ki + y)-%{ki)-\ 
T^{ki) = {^n(ki)\^n{h ^ fi)). (5) 

and subsequently the Chern number C = i^^^ ^xy{ki). 
Remarkably, the sequence of Hall plateaus is extremely 
sensitive to the values of the non-Abelian fluxes. In the 
Abelian regime <I>ck = = 0, we observe that the Hall 
conductivity follows the usual integer QHE cjxy = 
where the factor 2 is due to color-degeneracy. Con- 
versely, in the 7r-flux regime ($q, = = 7r/2) illustrated 
in Figj2jb, we obtain a completely different sequence of 
Hall plateaus where dxy = ^ ^) ^^o™^ = 0, 
as recently observed in graphene [4J. This sequence is 
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characterized by sudden changes of sign across the VHS 
situated at £^ = ±2, and by unusual double steps which 
can be traced back to the underlying low-energy relativis- 
tic excitations. As the gauge fluxes vary in the vicinity 
of the TT-flux point = 7r/2 + e and = 7r/2 — e), 
the system enters the non-Abelian regime and the Hall 
plateaus are modified (see Figj2]d). Indeed most of the 
degeneracies induced by the Dirac points are lifted and 
the anomalous double steps around = are progres- 
sively destroyed. However, a striking behavior occurs: as 
the non-Abelian fiuxes are varied, the two VHS originally 
situated at £^ = ±2 in the 7r-fiux point are split into four 



VHS 

red 



±2(1 + COS ^;3), E 



VHD 

green 



±2(1 + COS ^«), (6) 



as illustrated in Figj2]c for e = 0.1. Surprisingly enough, 
anomalous double steps in the plateau sequence reappear 
at higher energies outside the two red VHS, while the 
green VHS induce a sudden change of sign (see Figsj2]c- 
d). It is interesting to note that the anomalous behavior 
persists in the high-energy regime and that this effect can 
be probed by varying the parameter ^p. The tempera- 
ture required to observe these plateaus should be smaller 
than the spectral gaps, namely T ~ 10 nK. 

To identify the non-Abelian features in this QHE, we 
introduce the Abelian fiux ^ in the Dirac Hamiltonian (|3| 
by minimal coupling p ^ p ± ^(— and obtain 

= {g-cr^Ci ± g-cr~a^) ± {g^a^a^ ± g^cr~a) , (7) 

where cr+ = |xi)(X2|, c^" = |X2)(Xi| are color-fiip opera- 
tors, g± = {cy±Cx){Bo/2)^^'^ ^ and are bosonic chiral 
operators listed in ^4 . In the isotropic limit ^_ = 0, 
the Hamiltonian consists of an anti-Jaynes-Cummings 
term, a well-known interaction in quantum optics [25] 
that leads to the usual relativistic Landau levels (LL) 
recently observed in graphene [2]. Conversely, in the 
non-Abelian regime g- ^ 0, the Hamiltonian becomes a 
simultaneous combination of Jaynes-Cummings and anti- 
Jaynes-Cummings terms, producing a new type of Lan- 
dau levels. These novel LL are obtained by means of a 
Bogoliubov squeezing transformation S'(C) = e2 ^ ^ 

with C = — tanh~^(^_/^+), leading to the energy spec- 
trum 



^LLL=0, E^ = ±^(2B^CxCy)n, n = l,2... 
and corresponding eigenstates 
|LLL) = |x2)>S^(C)|vac), 

1^^) = ^IXi)^UC)|n - 1) ± Ux2)S^{C:)\n), 



(8) 



(9) 



V2I 



with |n) = (n!)~-^/^(a''')"'|vac) being the usual Fock 
states. Accordingly, the effect of non-Abelian fields is 
to squeeze the usual LL. In particular, the lowest Lan- 
dau level (LLL) is a zero-energy mode characterized by a 
colored squeezed vacuum, which is in clear contrast with 
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(c) Squeezed wavefunction 
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Figure 3: Vortex-like single-particle wavefunctions of the LLL 
0lll(^5 2/) for m = 4. (a),(b) Isotropic limit — Cy. (c),(d) 
Anisotropic regime Cy — 2cx. Note that distances are mea- 
sured in units of the magnetic length /b- 



its Abelian counterpart, the latter being simply the vac- 
uum. Besides, this LLL presents half the degeneracy of 
the remaining excited states n > 1 [26], and leads to the 
so-called anomalous half-integer QHE 



(10) 



where the filling factor v is defined as the integer part 
of [EY/2BQCxCy\^ and g is the Dirac points degeneracy. 
Let us stress that the non-Abelian fiuxes modify the Hall 
plateaus in a non-trivial manner as already emphasized 
through the numerical results. In particular, the Hall 



conductivity in Eq. (10) predicts the anomalous half- 
integer plateaus represented in Fig. [2|b), where the con- 
ical singularities are four-fold degenerate g = A. Con- 
versely, in the non-Abelian case shown in Fig. [2|d), the 
degeneracy is lifted to = 1, and thus the size of the 
steps is modified in accordance. 



As discussed above, the anomalous QHE is essentially 
a single-particle phenomenon that relies on the peculiar 
properties of the LLL. Additionally, further non-Abelian 
anomalies can also be found at the many-particle level, 
where an exotic Laughlin wavefunction [27j can be ob- 
tained by filling the single-particle vortex wavefunctions 



CLL(2^,y) = 



(11) 



Here Tx = Ib\/^^x/~c^j = ^b\/2c^7^7 describe the 
anisotropic extent of the wavefunction in units of the 
magnetic length Ib = \/T/Bo^ and m = 0, 1... represents 
the number of left-handed quanta [24 . Note how the 
loss of rotational invariance caused by the non-Abelian 
induced anisotropy Cx 7^ Cy^ leads to the squeezing of 
the vortex levels (Figs. [3]^a)-(d)). Filling these squeezed 
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degenerate states (11) according to Fermi statistics, we 



obtain the Laughhn wavefunction 



j<k 



(12) 

where u = cosh^, v = sinh^, f{u^v) = + v"^) 

and g{u^v) = \uv depend on the anisotropy through 
the squeezing parameter and Zjk = Zj — Zk repre- 
sents the complex two-fermion distance. In the Abelian 
limit C = 0? recovers the standard integer Laughlin 
wavefunction ^[z] = Ilj^/c /(^j' where 
f{zj^Zk) = Zj — Zk belongs to the space of holomorphic 
functions (Bargman-Fock space [28 ). Strikingly, in the 
non- Abelian scenario C 7^ 0? the wavefunction (12) does 



not belong to such space due to the interference between 
holomorphic f{z) and antiholomorphic f{z) components, 
and thus represents an instance of a non-chiral QHE. 
As shown below, this new anomaly modifies the clas- 
sical analogy with the one-component plasma (OCP), 
the building block that characterizes the peculiar prop- 
erties of quasiparticles in the fractional QHE [20^. The 
Laughlin state can be interpreted as the partition func- 
tion of a OCP (X Zc = JUidzjdzje-^^/^^ with 



kT 



1/2, a classical gas of particles interacting with a 



charged background through the potential 



C/e = -^l0g 

jk 



\uzjk-vzjk\ + \ J2j {f\zj\^ - giz] + z^)) . 

(13) 



The last term corresponds to the charged background 
jellium pj = —^^{^ + whereas the first describes 
a collection of positively charged particles q = I sur- 
rounded by a charge cloud Sp{z)^ with z = |z|e~^^, and 



Spi\z\,e) = 



tanhC (1 + tanh^C) cos 29 - 4tanhC 
\z\^ (1 +tanh\) - 4tanhCcos26'" 



(14) 



Notice how the surrounding charge cloud is absent 
Sp{z) = in the Abelian limit C = 0^ and we recover the 
usual OCP analogy. Conversely, for non- Abelian regimes, 
the collection of interacting positively charged particles 
becomes locally surrounded by an anisotropic charge 
cloud p = qS{z — Zj) -\-Sp{zj) with J d?z5p{z) = 0. In 
accordance, the paradigmatic plasma analogy is altered 
due to the squeezed nature of the LLL, a fact that may 
find profound consequences in the fractional QHE. 

We have shown that non- Abelian optical lattices offer 
an intriguing route to probe the striking properties 
of emerging Dirac fermions in anisotropic Minkowski 
space-times. In particular, the versatility offered by 
such experimental setups leads to the unique possibility 
of tuning the anisotropy of the underlying space-time, 
leading to remarkable effects such as non-chiral quantum 
Hall effects with several types of anomalies. 
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